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Abstract 

In this paper, we show that the Hamiltonian approach to loop quantum gravity 
has a fermion doubling problem. To obtain this result, we couple loop quantum grav¬ 
ity to a free massless scalar and a chiral fermion field, gauge fixing the many fingered 
time gauge invariance by interpreting the scalar field as a physical clock. We expand 
around a quantum gravity state based on a regular lattice and consider the limit where 
the bare cosmological constant is large but the fermonic excitations have energies low 
in Planck units. We then make the case for identifying the energy spectrum in this ap¬ 
proximation with that of a model of lattice fermion theory which is known to double. 
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1 Introduction 

Theories with a discrete notion of space have been a central theme of quantum gravity 
research since its birth. Whether assumed, as in Regge calculus or causal set models, 
used just as a regulator to be removed, as in causal dynamical triangulations, or derived, 
as in loop quantum gravity, in many different approaches to quantum gravity space or 
spacetime are represented by a discrete structure. However, there is a crucial problem 
with quantum field theories represented as degrees of freedom on discrete lattices, which 
is that it is impossible to represent chiral fermions by a local Hamiltonian. By a chiral 
fermion we mean one where the representation its left handed fermions live in is not the 
complex conjugate of the representation within which the right handed fermions trans¬ 
form. When we attempt to do so, we find that the spectrum of fermonic excitations dou¬ 
bles, with the right handed spectrum being now the mirror of the left handed one. 

We know this from examples, which are supported by a powerful no-go theorem 
proved by Nielsen and Ninomiya 1 for local quantum field theories on a regular lattice 
[1, 2, 3]. The result of this theorem is that for such models to be consistent, they must 
couple equally to the left and right handed sectors of the theory. As an immediate conse¬ 
quence, the number of left and right handed particles must be matched, and this result is 
referred to as fermion doubling. This theorem makes the construction of a lattice theory 
for the standard model quite difficult, since the weak interactions only couple to the left- 
handed sector of the the theory; if a lattice theory is to be compatible with the well tested 
experimental evidence of the standard model, it must avoid fermion doubling. 

1 For the interested reader we review the proof of this theorem in an appendix. 
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To date, most research into fermion doubling has been focused on standard quantum 
field theories. In this paper, we propose that fermion doubling is not unique to quan¬ 
tum field theories, rather we will provide a demonstration of fermion doubling in loop 
quantum gravity. 

One might expect that fermions double in loop quantum gravity just because space is 
discrete, but there is a deeper reason having to do with the absence of chiral anomalies. 
In a now classic series of papers [4, 5, 6, 7], Thiemann constructed a regularized form 
of the Hamiltonian constraint for quantum gravity-either the pure case or the case in 
which arbitrary collections of Yang-Mills fields, scalars and fermions are included. In all 
cases, Thiemann shows that the algebra of quantum constraints is first class and finite, 
hence, anomaly free. But consider the case in which a chiral gauge field is coupled to a 
multiplet of chiral fermions which have the feature that, when considered just as a chiral 
gauge theory, the dynamics would be inconsistent due to a non-vanishing chiral gauge 
anomaly. How can coupling to gravity render that inconsistent quantum theory finite 
and consistent? The only way this is possible is if fermions double, cancelling the chiral 
anomaly. 

More specifically, we work in the Hamiltonian formulation of loop quantum gravity, 
coupled to a scalar field and a chiral, left handed fermion field. The chiral fermion may 
also be coupled to a Yang-Mills gauge field. We follow [8, 9] and gauge fix the refoliation 
gauge symmetry by picking the scalar field as a time variable; this gives us a Hamilto¬ 
nian, U^ G . We expand the action of that Hamiltonian around a background state of the 
gravitational field which under coarse graining corresponds to a flat spatial metric. That 
background state is constructed from a superposition of spin network states with support 
on a regular lattice, T. We then define a Born-Oppenhiemer expansion around that state 
in addition to an expansion in the inverse of the bare cosmological constant. (These kinds 
of expansions were studied previously in [10, 11, 12] and [9, 13].) We then show that to 
leading order, the low energy spectrum of V, QG (for E <C E P i anck ) is equal to the spectrum 
of a lattice fermion theory on T. But we can use the Nielsen-Ninomiya theorem to show 
that the spectrum of that theory is doubled. Therefore, the fermonic excitations of the 
loop quantum gravity state also double. 

In a bit more detail, the basic strategy of our analysis is the following: 

1. Start with a lattice fermion theory based on a graph T embedded in E, a torus T d . 
Assume T is regular enough for a fourier expansion of the fermion excitations to 
exist, allowing there to be a continuous spectrum E = E(p). Verify the assumptions 
of the Nielsen-Ninomiya no-go theorem hold for the Hamiltonian H^ T . 

2. Define T-L®' G to be the (spatially diffeomorphism invariant, with a physical hamilto- 
nian coming from gauge fixing) subspace of the Hilbert space of LQG based on a 
graph T. 

3. Define H matter to be the Hamiltonian based on graph T that dominates in the Born- 
Oppenheimer approximation. 
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4. Use a degra vita ting map and dressing map to show that the spectra of and 

T must be identical, at least in the regime of E <C E P i anck . 

5. This establishes the no-go theorem for LQG. 

The moral of this strategy is to carry the problem of doubling in LQG to a standard 
lattice fermion theory, where the no-go theorem may be applied. Item 2 grants us a back¬ 
ground graph to work on. We use the Born-Oppenheimer approximation to concentrate 
on the fermionic degrees of freedom and expand around the gravitational pieces, we need 
the fermions to live on a background graph T to investigate doubling. Note our result is 
that there exists T such that fermion excitations of the corresponding background states 
double their spectra. We don't yet show that the fermion excitations around all states 
double, that would be a stronger result. 

2 Fermion Doubling in Loop Quantum Gravity 

In this section, we will work out the nontrivial items 2 to 4 of the program introduced in 
section 1. The main trick is to perform a Born-Oppenheimer like expansion [14] around 
a background state which is a functional of the gravitational degrees of freedom. We 
choose to work with a state with support on the spatial diffeomorphism class of a single 
fixed graph T, which we take to be a regular lattice (we do not expect our results to change 
if working with a state with support over a countable collection of graphs). By doing so 
we map the spectrum of low lying fermonic excitations of the quantum theory to that 
of a regular fermion theory on a fixed graph, which is the same T. We thus show the 
spectrum of this theory is equivalent to that of a lattice fermion theory which doubles. 
To demonstrate this equivalence, we will define degravitating and dressing maps, which 
carry the Hilbert spaces and Hamiltonians of the two theories into each other. In the low 
energy limit, the spectra of the two theories become equal. 

To reduce quantum gravity to lattice fermion theory, we shall work in the Hamiltonian 
framework of LQG, and study a model where quantum gravity is coupled to a scalar 
"clock" field T with conjugate momenta 7r, plus chiral fermion fields. 

We also extend the results to include a lattice Yang-Mills theory with compact gauge 
group G, in which case the fermions are extended to a multiplet of fermion fields in a 
representation r of G [15]. 

The fermion fields are represented by two component spinors ip Aa with conjugate 
momenta 7T 4 q , where A — 0,1 labels components of a Weyl spinor and a, which may 
sometimes be suppressed, labels the basis of the representation r 0 . For this theory, the 
Hamiltonian constraint, Diffeomorphism constraint, and Gauss constraint are a sum of 
terms 


3 



A, = Dr + Dl + D* + DZ m , G B b = G*” v + G* b , 


G — C*grav + Ct + + CyM , 

with 

Of = A? = G* AB = 7T (A ^ Bh (1) 

and 

C't = ^vr 2 + ^E?E M d a Td b T, D T a = nd a T, (2) 

= ^^«6(e°e 6 + 6 a 6 b ), £>™ = e afcc e & 6 c , (3) 

where all constraints with the superscript "grav" are the gravitational counterparts from 
Hamiltonian general relativity [16]. 

There is also the Gauss's law constraint of the Yang-Mills gauge group, G. 

G = V a e a + t(4) 

where e a and b" are the densitized electric and magnetic fields of the Yang-Mills field a a 
and here and above we have suppressed the indices for the fermion representation of the 
gauge group, r. 

To map a quantum gravity theory to a conventional quantum theory, we need to ex¬ 
tract a Hamiltonian from the Hamiltonian constraint. To do this we gauge fix the many fin¬ 
gered time refoliation invariance, before quantization. This gauge fixing is accomplished 
by taking the scalar field T to be a physical clock, so the spacial slices have constant 
T(x) = r. This turns the Hamiltonian constraint into a Schrodinger-like equation [8] 

ih^ = f d 3 xy/-2[C siav + ( 5 ) 

We will define W to be the regulated limit of the right-hand side 

w = j™ [ rf 3 x v t2[C“ + cy 5 ]. (6) 

The main work we have to do is to define W acting on spatial diffeomorphism invari¬ 
ant states. These are states which are functionals of diffeomorphism equivalence classes, 
which are labelled by the diffeomorphism class of the graph T, which is denoted by {T}. 
The state carries additional labels associated with the edges and nodes of T. Such states 
are of the form, 

*[T] = m = *[{r}] (7) 

The regulated operators involved in W are not diffeomorphism invariant, and so for finite 
values of the regularization parameters, the action has to be defined on a larger, kinemat- 
ical Hilbert space, 'H km . But, as discussed in detail in [17], the limit as the regulators are 
removed of the operators at issue give us an action purely on "H d d/ eo . 
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2.1 Classical Born-Oppenheimer Approximation 


We employ a state dependent, graph preserving regularization, so that the action of 
W on a state supported on a single graph T is another state supported on T. When acting 
on a state with support on a diffeomorphism equivalence class of a graph {T}, the action 
returns a state with support on that diffeomorphism equivalence class [18]. 

2.1 Classical Born-Oppenheimer Approximation 

In this section, we will expand the square root in the operator (12). To do this, we separate 
out the cosmological constant term 

dgrav = — det(T?)A + Cninst, (8) 

and perform an expansion around large A. Initially, this may look like a very poor ex¬ 
pansion. A is the infrared cutoff [19, 20] of the theory, and expansion around a large A 
suggests we are taking the infrared cutoff to be larger than the ultraviolet cutoff. How¬ 
ever, this is an expansion around the bare cosmological constant, which we already know 
to be large by the standard naturalness arguments. In addition, the particular value of 
A takes no role in our proof (see section 2.5), and we can add or subtract an arbitrary 
constant from it at will. 

The expansion of W for large A is now (with e as set of regularization parameters 
[8,15]) 


W = ! U o i d X \ Adet (^) - 2 [Qinst( a; ) + QXaH 


= lim / d 3 x 


e—>0 


\ 


Adet (E) 1 - 


2 C 


Ernst 


X) 


2 C^(x 


\ 


A det(-E) Adet(£) J 

= lim [ d 3 xV A v/det (E) - - 1 ^ X \ + 

A det (£')2 A det (£')2 


O 


i) 


= VXv-w eff , 

where the effective Hamiltonian W e ^ is 

yye// _ j- m 


& tx L d x 7m) ^ st(x) ++ 0 v ^ 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 


We can formulate the extension of the Nielsen-Ninomiya no-go theorem to LQG by 
establishing a map between the Hilbert spaces of two theories: LQG^, loop quantum 
gravity with chiral fermions (and possible Yang-Mills fields), and the lattice theory for 
fermions without gravity. 
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2.2 The Quantum Gravity Hilbert Space 


2.2 The Quantum Gravity Hilbert Space 

In this section, we shall review some key properties about the quantization. For simplicity, 
we shall fix the manifold E to be the torus T 3 . The Hilbert space of loop quantum gravity 
coupled to fermions and possible Yang-Mills fields will be denoted 

n QG *. (14) 

This is the spatially diffeomorphism invariant physical Hilbert space based on gauge fix¬ 
ing to a constant clock field T. On this Hilbert space, we will impose the spatial diffeomor¬ 
phism and G constraints. has a physical and spatially diffeomorphism invariant 

inner product. 

V,Q Q '* is decomposable in terms of diffeomorphism invariant classes of embeddings 
of a graph T into E 

n QG * = (15) 

Each component consists of extended spin network states, with labels correspond¬ 

ing to gravity, fermions, and Yang-Mills fields, diffeomorphic to T. One basis to the 
Hilbert space is the span of all spin-network states [21] of the form 

1.7 i, r, c, 0), (16) 

where j is the spin of an irreducible representation of the gravitational SU (2) associated 
to each edge of T, i is an 577(2) intertwiner associated to each node of T, while r and c are 
representations and intertwiners of the Yang-Mills gauge group G associated with edges 
and nodes respectively. Each node of the graph carries a basis state of the Hilbert space 
for fermions at a site, carrying a finite dimensional representation of 577(2) x G, spanned 
by the n and -0 operators. 

(Technical point, we impose equivalence under point wise smooth diffeomorphisms, 
so that there are no invariant diffeomorphism invariants characterizing classes of inter¬ 
sections with valence 5 or greater.) 

On , there is a physical Hamiltonian corresponding to the classical operator W 
(given by eq. (12)), obtained by gauge fixing the clock field to the T = constant gauge. 
We must now choose a regularization scheme. For simplicity we shall consider a graph 
preserving regularization, where W decomposes into separate actions on each graph sector 

W = J2W. r, (17) 

r 

where each Wr acts only on the corresponding H GC ' 1 '. We expect the results to be inde¬ 
pendent of this choice. We shall talk about this regularization choice in section 2.4 

Also defined on each are algebras of observables, 0® these are defined by 

loop and flux operators for SU( 2) and G, together with the fermion operators X, Y , etc, 
restricted to T. 
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2.3 Mapping Quantum Gravity to Lattice Gauge Theory 


Given a choice of graph T, there is also a lattice gauge theory involving the same 
fermion and Yang-Mills fields. The basis states are obtained by simply removing the 
gravitational degrees of freedom, 

\r,c,ip). (18) 

The lattice gauge theory Hilbert space 'H GT , contains an algebra of operators, 0 ( Y rr , an 
inner product, and a Hamiltonian H Y T . 

To summarize, the gravitational, gauge and fermion fields are all represented by a 
decorated spin network T, whose edged are labelled by (j, r ), representations of SU (2) x G, 
and whose nodes are labelled by intertwiners coming in one of the two types. In the next 
section, we will discuss the mapping between these two theories. 

2.3 Mapping Quantum Gravity to Lattice Gauge Theory 

In this section, we shall define maps between the two theories for each T, with the goal 
being to identify the spectrum of LQG based on a graph T with a lattice gauge theory. The 
simplest map to define is the degravitating map 

Gr : -»• U$ T , (19) 

which simply removes the labels corresponding to the gravitational fields 

Gr° \j,i,r,c,ip) = \r,c,ip) • (20) 

To go in the other direction, we define dressing maps 

Jr : (21) 

which takes a state of the lattice gauge theory and dresses it with gravitational fields; this 
is defined by a choice of amplitudes (pr[j, P r, c, ip] so that 

Jr o |r, c,ip) = Y^ <h\j, P r, c, ip] \j, i, r, c, ip) . (22) 

The relevant choice of these amplitudes will be clear in section 2.5. 

We require that 

Gr ° Jr = / (23) 

for the degravitaing map to be the "inverse" of the dressing map (Gr is onto, whereas Jr 
is into a subspace). 
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2.4 Regulating the Hamiltonian 


2.4 Regulating the Hamiltonian 

We need to define regularizations of the different terms in the Hamiltonian on the fixed 
graph subspaces (these are different than those previously mentioned since we are 
regulating on a fixed graph). 

Let us start with the fermion term in the time-gauge fixed Hamiltonian eq. (1) 

H*= [-C*= [ -n a A E aA {V a ^) B a . (24) 

J se 7s e 

Acting on a state <f> r in r we want to regulate this as 

H* H% g t=^Y. n (”)5 4 )I V(». i) B C V(n,a) a “<l(n + a) c0 , (25) 

^ (n,a) 

with n being the set of nodes and a being the corresponding adjacent vertices to each of 
these nodes. Both are summed over. It is the specification of these nodes and edges that 
make the regularization graph preserving. 

Here E(S)ab is the E\ B smeared over a two surface S with base point p defined by 

E(S) AB = j d 2 S( n \ n \E[S(a)Y CD U( lr ^)SU(r r \,^)S, (26) 

where is an arbitrary, non intersecting curve in S connecting the base point p with 

the point S(a). We pick the surfaces S(n, a) adopted to the graph T so that for each node 

n and adjacent edge (n, a), S(n, a) is a surface that crosses that edge once infinitesimally 
close to the node n base pointed at the intersection point. In addition, U(n, a) r f is the 
SU( 2 ) parallel transport over the edge (n, a), and V ( n, ajj 3 is the same for the Yang-Mills 
gauge group G. The specific action of the curves and surfaces are crucial to the definition 
of the regularization. 

We also have to define the fermion momentum operator. Let be the finite dimen¬ 
sional Hilbert space of fermions at the node n, which has a momentum operator 77 ', we 
have, when acting on a state in "Hj ? 1 , 

d 3 x s A (x)fl a (x) = J2 sA (n)nl (27) 

n 

The graph preserving regularization of the Yang-Mills part of the Hamiltonian works 
the same way. We do not give the details here but the important result is that the magnetic 
field squared term translates to a sum over minimal loops of the product of parallel trans¬ 
ports around those loops. A minimal loop is a c lose path on T that starts from and arrives 
to a base point, that is not made up of products of smaller loops. These are triangles in a 
simplical lattice and plaquettes in a cubic lattice. 
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2.5 Lattice based Born-Oppenheimer Approximation 


2.5 Lattice based Born-Oppenheimer Approximation 

We now apply the tools discussed earlier in this chapter to identify loop quantum grav¬ 
ity with lattice gauge theory As mentioned in the outline, we need to perform a Born- 
Oppenheimer approximation to expand around the gravitational degrees of freedom and 
work on one particular graph. To do this, we need to define the appropriate background 
state. Since T is fixed, we start by defining a class of states which have support only on T 

^r[A ,i,j] = (0, T|A, i,j) = Sr A f{i,j)- (28) 

We must define the background state \I/ 0 [T] in this form. We must now choose the particu¬ 
lar dependence of the state on the spins j, i. The dependence on the intertwiners i and j is 
chosen so that the state is an eigenvector of the Hamiltonian constraint with the smallest 
positive energy 

Mn jT^I°' r > = - B "“' |0 ' r >' (29) 

det(£) 2 

In addition, to replicate translation invariance, we will quantize the volume operator 

v/det(£)|0,r) = v. (30) 

The wave function which satisfies these properties will be referred to as To [T]. As is done 
in atomic physics [14], we use this choice of background state to perform the expansion 

T[r] = e^’To[r]y[r], (31) 

where the eigenvalue problem for T 0 fixes T, i,j, c is the coupling constant of the theory, 
and % is a wave function over the fermionic degrees of freedom |r, c, ip). We point out a 
subtle issue with the cosmological constant here: we may alter A and v in such a way as to 
preserve the following result. The value of v is of no physical importance to us, so instead 
of taking the IR cutoff to be larger than the UV cutoff, we may add appropriate powers 
of the UV cutoff to the cosmological constant without distrubing any of our results. To 
first order in the Born-Oppenheimer approximation, we neglect , and the resulting 
Schrodinger equation on a background state is 

lhd X [{T},r} = ^matter^-p}^ 

with 

H “' = +wx <0 ' r| L d3x vmi r> 

The constant E min may be subtracted without changing the results. At this stage, we 
perform the regularization in section 2.4. We evaluate the regulated operator on dif- 
feomorphism invariant states and take the limit as the regularization parameters are re¬ 
moved. As discussed in [18], the action is finite on these states and is within the Hilbert 


(32) 

(33) 


lim(0,T| f 

e_5>0 J s 
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2.5 Lattice based Born-Oppenheimer Approximation 


space associated with the graph T. The resulting LGT Hamiltonian is 

Hr GT = 77 ( n , + a) B g, ( 34 ) 

hl (n,a) 

where vQ & is tangent to the edge (n, a) at n and unit in the background metric, 

Qab = °a B °bAB, ( 35 ) 

as can be seen from the action of the Y a operator on a fixed background. We have now 
demonstrated one direction of the equivalence of the LQG and LGT theories, namely the 
degravitating one. We see the lattice spacing, a, comes out to be the Planck length. 

The same method applied to the regularization of the Yang-Mills contribution to the 
Hamiltonian constraint yields the lattice Yang-Mills action, with the lattice spacing com¬ 
ing out as the Planck length. 

We must now show the other direction, the dressing map, gives the same spectrum. 
Once this is complete, we will have shown these two theories are equivalent, hence LQG 
doubles. 

To perform the mapping, we need to define, for each graph subspace Tffr , the appro¬ 
priate degrees of freedom for the dressing map 

U ? G . ( 36 ) 

We insist the dressing state <f>p satisfies 

(4>°|£[S(n,a)]^|$°) 

and 

<*°|tf(n,a)/|<l>°} 

In addition, we also require 

«i-i*?>=4= 1 < 39 > 

e e u 

to ensure the volumes of cells are equal, as before. Notice we are setting a tensor density 
with weight one equal to a constant, thus this will only work in some frames. 

We note that the SU (2) Gauss's law constraint is not satisfied. 

Once we have defined the dressing and degravitating maps with the above choice of 
amplitudes, it is easy to show the spectra between the Born-Oppenheimer expanded LQG 
and the corresponding Lattice Gauge theory are identical 

<X| (*°r\HZ 9 rl*?> lx) = (X|^ GT |X), ( 40 ) 

with the lattice spacing taken to be l P[ . 


<,a°a B 


= 5b°, 


(37) 

(38) 
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This concludes the proof. In summary, we have performed a number of approxima¬ 
tions to transform a theory of loop quantum gravity on a background graph into a lattice 
gauge theory Upon reading this proof, one may worry about issues such as backreac- 
tion. However, it is critical to notice fermion doubling is a low energy phenomena. The 
left and right handed "limits" of a particle only make sense when the lattice spacing is 
taken to zero, in LQG the case is identical as the Planck length approaches zero. What 
makes the proof in LQG a nontrivial exercise is the notion of a background graph T: this 
is a purely gravitational phenomena and disappears when the Planck length is taken to 
zero. We must use the background graph for the question of fermion doubling to take any 
meaning. In addition, the inner product is not the most clearly defined object here, and 
the concept of momentum space is rather vague. These two points justify our expansions 
around a background graph and hint towards a mapping to a lattice gauge theory where 
we know how fermion doubling arises. 


3 Outlook 

In this paper, we have shown the spectrum of loop quantum gravity coupled to free 
fermions and Yang-Mills fields doubles. This provides an obvious contradiction with 
well known particle physics, as the weak interactions only couple to the left handed sec¬ 
tor of the standard model. In order to connect loop quantum gravity to well-established 
observations, we need to modify some of the assumptions put into this work. 

The main assumptions can be summarized as follows: 

1. Assume the bare cosmological constant A is large enough so that the square root in 
the Hamiltonian constraint can be expanded. 

2. The Born-Oppenheimer approximation is a reasonable assumption: The variation 
of gravitational fields does not effect the statement of fermion doubling. 

3. Assume the interactions between the fermions are local and Hermitian. 

4. Assume the momentum space is compact. 

5. Assert chiral invariance has the same definition for both lattice based theories and 
continuum theories, namely {D, 75 } = 0, with D the Dirac operator. 

6 . Assume the action is quadratic in the fermionic fields. 

We should emphasize that our result does not apply to spin foam models, which 
are path integral quantizations of loop quantum gravity. The reason is that the local 
Lorentz group is exactly realized and is represented by infinite dimensional representa¬ 
tions. Hence, the space of fermionic excitations is non-compact. 

Nor does our result rule out the possibility that there may be other backgrounds states 
of the quantum gravitational field, expansions around which lead to a chiral fermion 
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spectrum that does not double. Gambini and Pullin have presented an example in the 
case of a 1 + 1 dimensional reduction of loop quantum gravity [22], We note that their 
example involves a continuous superposition of background lattices, and will not apply 
to any superpositions of a countable number of (diffeomorphism invariant classes) of 
graphs. 

This highlights a missing element from loop quantum gravity research, which is the 
absence of a Hamiltonian for the gravitational sector, which is bounded from below [23]. 
Without this, we can construct semiclassical states which course grain to flat space, but 
we have no method to determine which is the best approximation to the ground state. 

Let us then consider how we might modify the Hamiltonian theory in a way that 
avoids fermion doubling. We believe the first four assumptions are reasonable. In par¬ 
ticular, fermion doubling is a low energy artifact of the theory, where it should be safe 
enough to ignore Planck-scale physics, and the cosmological constant can be added and 
subtracted at will. Instead, we suspect the problem is due to the choices related to the 
matter sector. Below, we will discuss several possible strategies to avoid fermion dou¬ 
bling. 

One of the requirements for the Nielsen-Ninomiya theorem to hold is that the fermion 
Hamiltonian must be quadratic in the fields. In standard QFT, this is a perfectly reason¬ 
able assumption to make as higher order terms are irrelevant; in 4D the coefficients of a 
higher order term in the fields has a positive mass dimension, hence they're irrelevant. 
However, in a lattice theory there is a natural cutoff provided, and in loop quantum grav¬ 
ity UV divergences are completely avoided, so in principle such terms can be added to 
the action at will. One possibility is that there may be a higher order action which does 
not double and is suitable for loop quantization. 

Another possibility to cure the problem would be to adapt known techniques from 
lattice gauge theory, such as the Ginsparg-Wilson technique, among others (these include 
Domain-Wall Fermions, Overlap Fermions, Twisted Mass QCD, Staggered Fermions, Wil¬ 
son Fermions, etc, see [24, 25, 26, 27, 28, 29, 30]). We could use the dressing map to define 
a Born-Oppenheimer expanded theory of quantum gravity which does not double, and 
use this to motivate a chiral regulation to the Hamiltonian constraint eq. (5). 

A third, and related, possibility is to look at perfect actions [31]. In this framework a 
lattice theory can exactly replicate the full set of gauge symmetries the continuum limit of 
the theory possesses. This is done through a process of coarse-graining the fields. How¬ 
ever, this comes at the price of locality. In addition, results can currently only be obtained 
in a perturbative framework, which is seemingly against the logic of LQG. 

Finally, we may hypothesize that chiral fermions are emergent states of the quantized 
gravitational field, such as the chiral braided states proposed in [32], In these states the 
chirality of the excitations is a consequence of the chirality of the braids themselves. 
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A Appendix: Review of fermion doubling in lattice mod¬ 
els 

A.l Review of Chiral Symmetry 

In this appendix we review the proof of the Nielsen-Ninomiya no-go theorem presented 
in [2]. Other, more rigorous proofs are listed [1, 3]. Before getting carried away with 
details, we shall take a small amount of time to discuss what fermion doubling physically 
means. 

In Quantum Field Theory courses [33], we are taught massless solutions to the dirac 
equation are eigenstates of the helicity operator 

h = p ■ S = pi ^ q ^ ^ , (41) 

with cr, being the Pauli Matrices. A particle with helicity h — +/- \ is referred to as a 
right/left handed particle. In other words, a right-handed particle's momentum vector is 
aligned with its spin (and anti-aligned for a left-handed particle). While this is perfectly 
consistent and clear in the normal continuum setting, in the discrete setting a particle 
must come with both helicities. To see this, we look at the fourier kernel e lpx ^' to notice 
a minimum length corresponds to a maximum momentum. Thus the momenta p and 
p + are identified, with a being the lattice spacing. We may add subtract off as many 
units of as we like, resulting in a positive momenta vector being identified with a 
negative one. Under this identification, the helicity changes sign and we see left and right 
handed particles must be equivalent. 

Another way to arrive at fermion doubling is to look at the dispersion relation in mo¬ 
mentum space. Since momentum space is bounded, any functions on it must be periodic. 
Thus the dispersion relation will cross the line of zero energy up and down an equal 
number of times. On each occasion, the continuum dispersion relation E = pc emerges 
in the low energy limit; however, it will come with varying sign and thus varying helic¬ 
ity. Since there are an equal number of zero energy crossings with positive and negative 
slopes, there are an equal number of emergent left and right handed fermions. As a re- 
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suit, fermion doubling occurs in all lattice theories with compact (closed and bounded) 
momentum space. This is the key feature we will look at in our proof. 

Helicity is not a meaningful observable if a particle is massive. This can be summa¬ 
rized quite simply, a massive particle does not travel at the speed of light, and hence we 
can consider an observer which overtakes the particle. In this frame, the direction of the 
momentum is flipped, and the helicity changes sign. For massive particles, it is conve¬ 
nient to think of a more abstract concept, known as chirality. The chirality of a particle is 
determined from whether the particle transforms under the left or right-handed part of 
the Poincare group. From here, we may work with chirally invariant actions, where the 
action has a symmetry under 


= e*^L, fn = e*N/; R , (42) 

where we can define the left and right handed pieces using the projectors 
Pr = While chirality is a reasonably defined concept for massive particles, it is easy 
to show mass terms break chiral symmetry. To do this, we can write a mass term using 
the left and right-handed parts 

rmjnp = m(-0z,'0L + i’L'^R + f’R'f’L + t/Wr)- (43) 

We see the cross terms break chiral invariance, and a chirally invariant theory must be 
massless. For a massless theory, chirality is conserved for a massless particle and is equal 
to its helicity. 

A.2 Nielsen Ninomiya no-go Theorem 

We are now prepared to formulate the no-go theorem: 

Theorem 1. Suppose we are given a lattice theory of free fermions governed by an action quadratic 
in the fields 


S = 


dt ^ ^ t^j I dt £ 


(44) 


x,y 


with f>(x, t ) an N component spinor with discrete position label x. Then fermion doubling is in¬ 
evitable if the following conditions are met: 

• The underlying lattice has a well defined momentum space which is compact. 

• The interaction H is hermitian and local, in the sense that its momentum space continuation is 
continuous. 

• The charges Q are conserved, locally defined as a sum of charge densities Q = = 

f(x) r f{x), and is quantized. 


14 




A.2 Nielsen Ninomiya no-go Theorem 


The proof below can be generalized to a more general kinetic term 'ip(y,t)T(x — 
y)ip(x,t) so long as the term T doesn't vanish anywhere on the lattice, as can be seen 
from the equations of motion. 

i^2 T (x- y)d t 'ip(y, t) = H(x - y)ip{y). (45) 

y y 

Inverting the T matrix, the equations of motion indicate the action is equivalent to one 
governed by the standard kinetic term along with the effective Hamiltonian H/T; if T is 
singular or nonlocal the no-go theorem does not apply Because of this it is often conve¬ 
nient to speak in terms of the Dirac operator D, defined by 

S = a 4 ifiDip, (46) 

X 

where now we add an additional assumption to the theorem: D is required to be invert¬ 
ible. 

In addition, we may couple to the fermions a fixed gauge field (such as a Yang-Mills 
field) without spoiling the proof; such fields simply act as a background when working 
in the low-energy regime. 

To characterize fermion doubling, we will look at the dispersion relation in momen¬ 
tum space, determined from the eigenvalue problem 

H(p)A(p) = w*(pM(p), (47) 

here the index i is not to be confused with the components of ip, but rather a label of 
eigenvectors which is not summed on. We can choose this label such that the eigenvalues 
are increasing (recall they are real by the Hermiticity of H) 

ui(p) < ca 2 (p) < • • • < u N (p). (48) 

The relevant feature will be captured by the level crossings 

ip deg) i+l{jp deg) i (49) 

as this corresponds to the zeroes of the energy difference in our simple picture. We will 
now show that the level crossings correspond to left or right handed Weyl particles. 

As an aside, notice that level crossings between 2 levels are generic in 3+1 dimensions, 
but 3 level crossings are not. To see this, consider a general 3x3 Hamiltonian spanned by 
the Gell-Mann matrices A, 

8 

H^\p) = A(p)l + J2 B i(p)^- (50) 

i= 1 

For a three level degeneracy the coefficients Bi(p ) must all vanish. This requirement gives 
8 equations for only 3 quantities, which generically has no solutions. However, two level 


15 




A.2 Nielsen Ninomiya no-go Theorem 


crossings are generic in 3+1 dimensions as the relevant 2x2 piece of the Hamiltonian is 
spanned by three Pauli matrices. 

Near a 2-level crossing at the degeneracy point p deg , we may expand the relevant piece 
of the Hamiltonian in a Taylor Series 

# (2) (Pdeg + Sp) = o;,:(pdeg)l + Sp ■ A + Sp k V^a a + 0(Sp 2 ) (51) 

where a and k run from 1 to 3. To simplify the form of this Hamiltonian near p deg , we will 
shift the momentum by 


P 0 = H -+ P 0 - Ui(p deg ) - Sp ■ A (52) 

5 Pk -+ S Pk ± S Pa V k a , (53) 

where the ± sign depends on the determinant of V. The shifted Hamiltonian becomes 

H = p-a, (54) 

and the corresponding eigenvalue problem is 

p . aU(p) = ±p 0 U(p ), (55) 

where U(p) is the corresponding wave function. The sign of the determinant of V is thus 
identified with the helicity of a Weyl particle. Thus each 2 level-crossing p deg is identified 
with a Weyl fermion in the low energy limit. It is important to note the helicity depends on 
the degeneracy point p deg and so it now remains to identify which types of Weyl fermions 
emerge in the low energy limit. 

To set about this task we will look at curves in the 3-D dispersion relation space (which 
is embedded in the full 3+1D u-p space) defined by 

{(P,w)| (a\uji(p)) = 0}, (56) 


where the bracket is 


(a\ui(p)) = aiV’P + a 2 'ip2 ) 4-4 (57) 

with | a) being any constant N-vector. In 3 dimensions, this will turn into 1 complex equa¬ 
tion for 3 quantities p, implying the set forms a curve. These curves are of special impor¬ 
tance because they pass through all of the degeneracy points p dcg . The reason for this is 
that at a degeneracy point, we may redefine the energy eigenstates through any superpo¬ 
sition of the form 

\u>i) = a\wi) +/3\u i+ i). (58) 

Thus we may choose a and f3 so that ( a\ui(p deg )) = 0. All such curves must be closed 
since the Brillouin zone is compact. Thus if the curves have an orientation, they must pass 
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equally many times up and down through the degeneracy points, proving the theorem. 
In the paragraphs below we will set up this orientation. 

To set about this task, we will look at the phase of (a\ui) on small circles wrapping 
around the curve near each degeneracy point. Suppose we set this small circle of radius 
R a distance d away from a degeneracy point atp 2 = 0 along a curve passing through the 
positive p z direction. The Weyl equation near the degeneracy point is given by eq. (55), 
with the sign being identified with helicity. The eigenvectors to this equation along the 
mentioned circle S 1 = {6 G R\(p x ,p y ,p z ) = (R cos 6, R sin 9, d)} in the limit R <C d are 

U\ = ( B_ e ie V wi = ±(d + R), U 2 = ( _R_ e -ie ^ , ^2 = =F (d + R). (59) 

V 2 dr J V 2 a J 

In the limit R —> 0, we have U\, 2 = ^ J ^, implying the vector (a| along the curve is (0,1). 
We obtain the following for the phase rotation along the curves for both levels 

(alt/,} = ±i-e* <a|t7 2 ) = (60) 

2 p z 2p z 

eq. (60) may now be used for orientation assignment. There are then two types of curves 
for each helicity, one crossing from p 0 < 0,p z < 0 to p 0 > 0,p z > 0 , and one crossing from 
p 0 < 0,p z > 0 to p 0 > 0, p z < 0. Using eq. (60), we see the first curve is oriented in the 
positive (negative) p z direction for a right (left) handed degeneracy point and the second is 
oriented in the negative (positive) p z direction for a right (left) handed degeneracy points. 
Since the curve must be closed and passes equally many times up and down through 
degeneracy points, by the orientation assignment we see the number of left and right 
handed particles are matched. □ 
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